Influence of gauge fluctuations on fermion pairing order 

parameter 



Yu-Liang Liu 

Center for Advanced Study, Tsinghua University, Beijing 100084, People's Republic of China 

Abstract 



Using a prototype model, we study the influence of gauge fluctuations on 
fermion pairing order parameter which has the gauge symmetry, and demon- 
strate that the gauge fluctuations can destroy the long range order of the 
fermion pairing order parameter, and make it only have short range correla- 
tion. If this parameter is a superconducting order parameter, we show that 
the Meissner effect of the system keeps intact, and the system is in the su- 
perconducting state even though the long range order of the superconducting 
order parameter is destroied by the gauge fluctuations. Our calculations sup- 
port that the pseudo-gap region of the high Tc cuprate superconductivity is 
a spin pseudo-gap region rather than an electron pre-paired region. 
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After the discovery of the high Tc cuprate superconductor |T| , it has been soon reached a 
common consensus that the high Tc cuprate superconductors are a two-dimensional strongly 
correlated system, and low energy behavior is determined by their copper-oxide plane (s) 
0-0. However, so far a complete theory representing this system does not appear, and 
there are still a lot of controversies and unsolved questions in this field. This strongly cor- 
related electron system has novel physical properties of normal and superconducting states. 
In normal state, the system shows non- Fermi liquid behavior and has strong anti- 

ferromagnetic fluctuations in low doping and optimal doping regions. Moreover, in the low 
doping region there is a pseudo-gap region, and stripe phase fluctuations may also appear. 
These phenomena cannot be clearly explained by an unified theory. In superconducting 
state, the system has a d-wave symmetry superconducting order parameter, but it is not 
clear which mechanism can produce this order parameter. There is another puzzle question 
that single electron excitation spectrum keeps intact as the system going to the supercon- 
ducting state from the pseudo-gap region PHTT|1 . This phenomenon seems to mean that the 
pseudo-gap region is an electron pre-paired region [|T2] . 



Usual perturbation methods are hard to treat strongly correlated systems, because in 
which there is not a controllable small parameter, and electrons have strong correlation. 
For enough strong electron correlation, there are not well-defined quasi-particles near the 
Fermi surface, and the systems show non-Fermi liquid behavior. On the other hand, these 
perturbation methods cannot directly describe the electron correlation. However, the elec- 
tron correlation is a key parameter hidden in the strongly correlated systems. The eigen- 



functional bosonization method |13| is a good candidate for treating the strongly correlated 
systems, because in this method there naturally appears a phase field, where its imaginary 
part represents the electron correlation. 

In this paper, we introduce a prototype model, which represents a two-dimensional spin- 
1/2 fermion system with transverse gauge fields, to study the influence of the gauge fluctu- 
ations on the fermion pairing order (spin pseudo-gap) parameter, here we assume that the 
fermion pairing order parameter does not violate the gauge symmetry, and these transverse 



gauge fields originate from some other strong fermion interactions, such as the Hubbard 
model with large on-site electron Coulomb interaction, or the electron single occupied con- 
straint of the t- J model. In the absence of the gauge fields, if the fermion pairing parameter 
is a constant (long range order), we demonstrate that after turning on the gauge fields, 
the correlation function and the average value of the fermion pairing order parameter go to 
zero with anomalous exponential asymptotic behavior. The transverse gauge fields provide 
a strong phase fluctuation factor to the fermion pairing parameter by phase fields. If these 
fermions are electrons, we demonstrate that the Meissner effect of the system keeps intact, 
and the system is still in the superconducting state, even though the long range order of the 
superconducting order parameter is destroied by the transverse gauge fields. This result is 
similar to that of a two-dimensional superconductor with long range Coulomb interaction, 
where the long range order of the superconducting parameter is destroied by the long range 
Coulomb interaction, but the Meissner effect of the system leaves intact |14| . However, there 



is a little difference between these two cases, for the long range Coulomb interaction, the 
correlation function of the electron pairing order parameter shows power-law asymptotic 
behavior. 

We consider a prototype model represented by the Hamiltonian with transverse gauge 
fields A(x) (V- A(x) = 0), 

H = J tfx^i*)^-™ - gAf - ^{x) 

+ V4(z)[^HV + gA) 2 - n]fo(x) (1) 
+ A(z)^i(z)V> 2 0e) + A*(z)^(ar)^(a:)} 

where the fermion field ipi{x) represents the fermions with spin-up, the fermion field ip2{x) 
represents the fermions with spin-down, and A(x) = — u < tpl(x)ipl(x) > is fermion pairing 
order parameter (spin pseudo-gap) induced by other interactions, where u is the interaction 
strength. It is noted that the fermion pairing terms do not violate the gauge symmetry 
of the transverse gauge fields. At A(x) = 0, for simplifying our following calculations, we 
simply take this parameter as a real constant, A(x)|a=o = A . The system described by the 



Hamiltonian (|]) is a strongly correlated system, as A(x) =0 it shows the non- Fermi liquid 



behavior P!5|,IB|. 

With the Hamiltonian ([]]), the action of the system can be written as a simple form, 
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where ^(x, t) = (ip\(x, t), tp 2 (x, t)), M = M + (f)(x, t) and 
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We have separated the propagator operator M(x, t) into two parts, one is a free term, and 
another one depends on the transverse gauge fields A(x,t). In this form, we first exactly 
solve the eigen-equation of the operator Mq, then we determine the differential equations of 
phase fields induced by the transverse gauge fields. 

Integrating out the fermion fields tpi(x,t) and ip2{x,t), we obtain the effective action 
(omitting constant term), 



S[A] 



d£ / dtd 2 xtr (ci>(x,t)G(x,t;x',l/, [£</>]) 



i f^t 

x'— >x 



(3) 



where the Green's functional G(x, t; x', t', [0]) can be represented by eigen-functionals of the 
propagator operator M(x,t), 



(4) 



The eigen-functional equation of the operator M(x, t) reads, 



M(x,t)^l(x,t, [0]) = ^j[0]^(x,t 



(OuhTrW, 



(5) 



where i — 1, 2 and the eigenvalues are E^' 2 ^] = wi-Efc+Ej. 1 '^ [</>], where = \Je 2 {k) + A 2 ,, 
e(Jfe) = fc 2 /(2m) -/i, and S^[0] = / X d£/ dtd 2 x¥£l*(x, t, [<j>])i(x,t)V$,(x,t, [(f)]) is a smooth 



(1,2) r 



function, and independent of u. Due to the non-zero fermion pairing parameter Ao, the 
energy spectrum of the fermions splits into two branches, one represents unoccupied states 
(conducting band), and another one represents occupied states (valence band). The fermion 



il 2) 

pairs are in the top of the valence band. The eigen-functionals *& k(J (x, t 
occupied and unoccupied states, respectively, and they can be written as 
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where u k = y(l + e(k)/E k )/2 and v k = y(l — e(k)/E k )/2 are coherence factors, A k is a 
normalization constant, and T and L are the time and space length scales of the system, 
respectively. At A(x,t) = 0, the eigen-functionals in (|[) are the exact solutions of the 
operator M , where A k — 1. The fermion pairing parameter now can be approximately 
written as near the Fermi surface (k ~ kp), 



A(x,t) = Aoe^^+qtCM) 



(7) 



where A is determined by the self-consistent equation A = 272 J2 k \A k \ 2 A /E k . The phase 



fields Q k (x,t) and Q k (x,t) satisfy Eikonal-type equations [fL3|l . If we only keep linear terms, 
we can obtain that, 



k • V V 2 a 
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Under this approximation, we have the same effective action as that obtained by usual 
random-phase approximation (RPA) (see below). By solving these equations, we have the 
relations between the phase fields and the transverse gauge fields, 
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where D ± = (O ± k • q/m) 2 - g 4 /(2m) 2 t A g 2 /m, Q fc ( 9 , fi) = Q%(q, fi) + Qj[(g, fi), and 
Qk(q, ^) = Qkili ^) + The imaginary parts of the phase fields do not contribute 

to the effective action, but they represent the fermion correlation induced by the transverse 
gauge fields. Only the real parts of the phase fields contribute to the effective action of the 
system. 

With equations (§]) and @, we have the effective action, 

s l A ] = ^E%(9.^H.-^ife n ) ( 10 ) 

where Iltffofi) = («r(A )fi/g + x(Ao)<? 2 )(% - <Mi/<? 2 ), and r(A ) = t(Ao) + 2mA (7 + 
7(A ))/g 2 . The parameters 7(A ) and x(Ao) are the smooth functions of the gap A . At 
Ao = 0, we have 7(0) = 7 and x(0) = X, an d have the same effective action of the transverse 
gauge fields A(x, t) as that of Refs. |||7[], where it was obtained by usual RPA method. We 
would like to point out that in the eigen-functional bosonization method it is naturally to 
introduce the phase field, which is a key parameter hidden in strongly correlated systems, 
because its imaginary part represents the fermion correlation. However, it is absent in 
previous perturbation methods. With the phase field, we can easily study the influence of 
fermion interactions on the low energy behavior of the system, and calculate a variety of 
correlation functions. 

At A = 0, the system shows non- Fermi liquid behavior produced by the strong fluctua- 



tions of the transverse gauge fields ||i5l , |l6f . It can be easily shown by present representation. 
With equation @ and taking functional average of the transverse gauge fields A(x,t), we 
can obtain single fermion Green's function, 

Gi(x — x ,t — t') = i < Gn(x, t] x , t', [(f)]) >a 

= J_ ^(-e(A;))e ik ' (x_x ' )_: ' e(A '' )( * _ *' ) e Pfc/(a: " a; ''*"*' ) (11) 

k 

where Pi (x, t) comes from the contribution of the imaginary part of the phase field Qk{x, t), 
and the contribution from the real part of the phase field Qk(x,t) can be neglected. Near 
the Fermi surface k ~ kp, we have the relation, 
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Pji r {0,t) ~ (12) 

where a = j^i^ix/j) 1 J dx(l — cos(x))x~^ 3 . This result is basical same as that of Ref. 
15| , p76[1 . Due to the time dependence of the phase factor P£ (0, t), the single fermion Green's 
function shows a singular low energy dependence which violates the quasi-particle excitation 
assumption of the Landau Fermi liquid. Therefore, the system shows non-Fermi liquid 
behavior. 

We now study the influence of gauge fluctuations on the fermion pairing order parameter. 
From equation (|7]), we see that the gauge fields A(x, t) provide a phase fluctuation factor to 
the fermion pairing parameter A(x, t) by the phase fields Qk{x,t) and Qk{x,t). Therefore, 
the gauge fields strongly affect the fermion pairing parameter. By taking functional average 
of the transverse gauge fields, we have the following relations, 

< A(x,t)A*(x,t) > A ~ Aq, < A(x,t)A*(x',t) > A ~ Age- 6 !*-*', 

< A(x,t)A*(x,t') > A ~ e - c| *-*' |1/5 , < A(x,t) > A ~ e- d/uJ o\ (13) 

where b ~ A 1/5 , c ~ 2 1 / 5 i ^[x(Ao)/(2mAo(7 + 7(A )))] 1/5 , 2d ~ c is a constant, and u is 
a characteristic energy scale of the transverse gauge fields, in general it is very small (ujq ~ 
1/L). These relations are our central results, we now explain them. The first expression 
means that the real parts of the phase fields Q/.(x,t) and Qk(x,t) have little influence on 
the spin pseudo-gap, and the last three relations mean that the imaginary parts of the phase 
fields strongly suppress the correlation of the fermion pairing parameter A(x, t), and make 
the long range order of the fermion pairing parameter be a short range order. Therefore, the 
fluctuations of the transverse gauge fields can destroy the long range order of the fermion 
pairing parameter, and make it only have short range correlation. 

If these fermions are electrons, the pairing parameter A(x,t) is a superconducting order 
parameter. At A(x,t) = 0, the Hamiltonian (fj) represents a two-dimensional superconduc- 
tor with superconducting gap Ao- The transverse gauge fluctuations make the system have 
a short range pairing parameter A (re, t), but they do not destroy the superconducting state 
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of the system. This can be justified by calculating the Meissner effect of the system after 
turning on an external magnetic field B(x,t) = V x A e (x,t). The current functional reads, 

J(x,t,[A}) = -ij{(x,t)— ifii(x,t) -4(x,t)— Mx,t) (14) 

m m 



- — [(A e (z, t) + gA(x, t))ip\(x, t)ipi(x, t) + (A e (x, t) - gA(x, t))ipl(x, t)if> 2 (x, t)} 
Taking functional average of the transverse gauge fields, we have the relation, 



J(x,t) = < J{x,t, [A]) > A 
fa{A ) p(A o y 



m? m 



A e (x,t) (15) 



where p(A ) = -(1/L 2 ) Z k 9(-e(k))e(k)/E k , and a(A ) = (1/L 2 ) E fc ^ ^gg^ l k ^ k - 
At A = 0, we have 3(x,t) = 0. This equation shows that even though the transverse 
gauge fields destroy the long range superconducting order, the Meissner effect of the system 
keeps intact, i.e., the system still is in the superconducting state. This result is consistent 
with that of equation ([13]). The Meissner effect is mainly determined by the coherence 
factors Uk and Vk of the superconducting state, and is insensitive to the phase factor of 
the electron pairing order parameter. However, the fluctuations of the transverse gauge 
fields can destroy the long range order of the electron pairing order parameter. One must 
remember that these gauge fields are not electromagnetic gauge fields, they originate from 
some strong electron interactions, and the superconducting state of the system does not 
violate this gauge symmetry. 

Even though all our calculations are based on the s-wave symmetry of the pairing pa- 
rameter, they are also qualitatively valid for the d-wave symmetry of the pairing param- 
eter. The above results may provide some useful informations for studying the high Tc 
cuprate superconductors, where their low energy physical properties are determined by the 
two-dimensional strongly correlated electron system in the copper-oxide plane(s). It is well- 
known that this strongly correlated system can be well described by the t-J model [Q or 
the Hubbard model with large on-site Coulomb interaction 0. Electron single occupation 
constraint of the t-J model in fact represents the strong electron correlation, which can 
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induce the gauge interaction among the fermions and holons in the slave boson/fermion 



representation of the electron operator [|1y]-|T9||, or the gauge interaction of the electrons 
due to the Hilbert space of the system is strongly suppressed by this constraint. We now 
only focus on the pseudo-gap region in the low doping region. If this region is an electron 
pre-paired region, and the change of electron pairing order parameter from long range order 
to short range order is induced by the gauge fluctuations, our calculations show that un- 
der an external magnetic field, the Meissner effect of the system will be observed, and the 
system is in the superconducting state, which contradicts with present experimental data. 
In fact, the system is in a normal state, and there is not the Meissner effect in this region. 
However, our calculations support that the pseudo-gap region is a spin pseudo-gap region, 
and the fermions do not have electric charges. The transport property of the normal state is 
determined by the holons (carrying electric charges). Due to strong gauge fluctuations, the 
system still shows non-Fermi liquid behavior even though there is the spin pseudo-gap in 
spin excitation spectrum. This is qualitatively consistent with present experimental obser- 
vations. How does the system go into the superconducting state? There may have two ways, 
one is usual Bose-Einstein condensation of the holons [0,171, however, in this case there is 



too high onset temperature of superfluidity. Another one is that in the low temperature 
limit the gauge fluctuations become enough strong to re-combine a fermion and a holon 
into an electron, and the fermion pairing order parameter becomes the superconducting 
order parameter. In this case, the superconducting order parameter has the same symme- 
try (d-wave symmetry) and the same modulus with the fermion pairing order parameter, 
which is qualitatively consistent with the angle-resolved photoemission experiments PHTT 



where single electron excitation spectrum nearly has no change as the system going to the 
superconducting state from the spin pseudo-gap region. 

In summary, with the prototype model, we have studied the influence of the transverse 
gauge fields on the fermion pairing order parameter which keeps the gauge symmetry, and 
demonstrated that the gauge fluctuations destroy the long range order of the fermion pairing 
order parameter, and make it only have the short range correlation. If this order parameter 



is the superconducting order parameter, we demonstrated that the Meissner effect of the 
system keeps intact, and the system is in the superconducting state even though the long 
range order of the superconducting order parameter is destroied by the gauge fluctuations. 
Therefore, our calculations support that the pseudo-gap region in the low doping region of 
the high Tc cuprate superconductivity is the spin pseudo-gap region rather than the electron 
pre-paired region, because the high Tc cuprate superconductivity is a strongly correlated 
electron system, where the Hilbert space of the electron states is strongly suppressed, the 
gauge fluctuations are existing. 
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